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The magnetic field dependence of the electronic specific 
heat C(H) in the s-wave superconductor NbSe2 shows cur- 
vature at low fields, resembling the near yH term in C(H) 
which has been reported in high-T c superconductors and at- 
tributed to a d x 2_ !/ 2-wave pairing state. In NbSe2 we find that 
the low-field behaviour in C(H) is described quantitatively by 
the expansion of vortex cores and the field dependence of the 
magnetic induction B above H a \. The associated change in 
the density of quasiparticle states localized in the vortex cores 
provides a simple explanation for the "low-field anomaly" in 
C(H) observed in s-wave superconductors. 
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In the vortex state of a type-II superconductor, an ap- 
plied magnetic field penetrates the bulk in the form of 
quantized flux lines. The predictions of this state by 
Abrikosov H and of bound quasiparticle (QP) states in 
the normal vortex cores of a conventional superconduc- 
tor by Caroli et al. Q remain two great achievements in 
superconductivity theory. Although flux decoration ex- 
periments readily confirmed the existence of the vortex 
state, the electronic structure of the cores predicted by 
Caroli et al. was not established until the advent of scan- 
ning tunneling microscopy (STM). In particular, Hess et 
al. |3|] observed tunneling spectra in the vortex cores of 
NbSe2 consistent with localized QP states — thus seem- 
ingly completing the picture of the vortex state. 

Interest in the vortex state was renewed with the dis- 
covery of the high-T c superconductors (HTSCs). The 
general consensus is that the charge carriers in the 
HTSCs form pairs whose wavefunction (or order param- 
eter) has a dominant d 2 .2_ !/ 2-wave symmetry, rather than 
the s-wave symmmetry characteristic of low-T c conven- 
tional superconductors. Much theoretical work has fo- 
cussed on incorporating d x 2_ y 2-wave symmetry into a 
model for the vortex state of the HTSCs. Despite these 
efforts, experiments have yet to unambiguously confirm 
these predictions. 

Measurements of the electronic specific heat C (T, H) 
are one way of probing the QP excitation spectrum in the 
vortex state. In an s-wave superconductor, where there 
is an isotropic energy gap at the Fermi surface, there 
is a contribution to C'(H) which is proportional to the 
QP density of states (DOS) localized in the vortex cores. 
Since the density of vortices increases linearly as a func- 
tion of magnetic field, this term is expected to be propor- 



tional to H Q] . On the other hand, in the vortex state of 
a c? a .2_j / 2-wave superconductor, Volovik ]jj predicted that 
the DOS varies as yH. This weaker field dependence of 
the DOS is mainly due to delocalized QPs which leak 
outside of the vortex cores along the nodal directions 
of the order parameter. Within this model, bound core 
states are at most a minor correction to the total DOS. 
Experiments performed on the HTSCs, YBa 2 Cu30y_,5 
(YBCO) §§ and La2-*Sr x Cu04 (LSCO) §,§, showed 
that there is a contribution to C(H) which is approx- 
imately proportional to y/~H — consistent with Volovik's 
prediction. Although this was one of the key early ex- 
periments providing evidence for d x 2_ y 2-wa,ve symmetry 
in the HTSCs, the interpretation of such measurements 
has been plagued by puzzling reports of similar curva- 
ture in C(H) at low magnetic fields in s-wave super- 
conductors, such as NbSe 2 (Ref. |o|), V 3 Si (Ref. Q) 
and CeRu2 (Ref. jl2)) — and in other unconventional su- 
perconductors, like the heavy fermion superconductor 
UPt 3 (Ref. the organic superconductor (BEDT- 

TFF) 2 Cu[N(CN) 2 ]Br (Ref. (TJ) and the borocarbide su- 
perconductor LuNi 2 B 2 C (Ref. fllSfl ), Ramirez sug- 
gested that this behavior at low fields might be a general 
feature of all superconductors in the vortex state, inde- 
pendent of the order parameter symmetry, but somehow 
related to the strength of the vortex- vortex interactions. 
Clearly the \H dependence of the specific heat in the 
HTSCs cannot be attributed to nodes in the energy gap 
function without a satisfactory explanation for similar 
behavior in fully gapped superconductors. 

A long-standing general belief is that the supercon- 
ducting coherence length £ is independent of H . How- 
ever, muon spin rotation (//SR) |T(| and STM jlTj mea- 
surements have shown that the radius ro of the vortex 
cores in NbSe 2 expand at low fields. In the /iSR experi- 
ment, the field dependence of ro could be fit to a function 
dependent only on the intervortex spacing. These /zSR 
and STM results have been met with some skepticism, 
however, since the change in ro implies that the super- 
conducting coherence length £ also varies with field. An 
obvious question raised by these observations is: Could 
the low-field anomaly in the specific heat be due to the 
expansion of the vortex cores? 

To address this question we have carried out detailed 
measurements of the specific heat C(T, H) in a 40 mg 
crystal of NbSe 2 from the same batch used in the /zSR 
experiment of Ref. . The crystal had a superconduct- 
ing transition temperature T c = 7.0(1) K, as determined 
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previously [ jl6[ by resistivity and susceptibility measure- 
ments, and here from measurements of the specific heat in 
zero field (see Fig. 1 ). From magnetization and specific 
heat measurements, the upper and lower critical fields 
for this crystal were found to be H C 2 = 2.90(3) T and 
H cl = 0.025(2) T at T = 2.3 K, and H c2 = 1.75(2) T and 
H c i =0.015(1) T at T = 4.2 K. The specific heat was mea- 
sured using a thermal relaxation calorimeter |l8f| with the 
magnetic field applied normal to the NbSe2 layer direc- 
tion. The data have been corrected for the small field- 
dependent background and addenda contributions. 

Figure 2 shows the field dependence of C(T, H) /T 
in NbSe 2 at T = 2.3 K and below H c2 . The down- 
ward curvature at low fields is similar to that reported 
m Ref. 0. Measurements were taken four different 
ways: (i) The single crystal was zero-field cooled (ZFC) 
to T = 2.3 K, and measurements taken for increasing field 
up to iJ = 2.5 T, (it) then taken for decreasing field down 
to H = T. (Hi) Starting with the crystal at T — 2.3 K and 
H = 4 T (i.e. in the normal state), the field was ramped 
down to 1 T and measurements taken for decreasing field 
down to H — T. (iv) Measurements were also made after 
field-cooling (FC) the crystal to T = 2.3 K, in different 
magnetic fields. No significant hysteresis is found in the 
specific heat measurements — consistent with the near re- 
versibility of the magnetization (Fig. 2, inset). 

In a conventional s-wave superconductor, the specific 
heat in the vortex state (H>H c i) is greater than in the 
Meissner state (H <H c i), due to a contribution from the 
localized QPs in the vortex cores. A precise calculation, 
from the Bogoliubov equations assuming noninteracting 
vortices, gives the density of states N(E) per unit volume 
associated with the bound excitations as B 



N(E) = 7V(0Ko 
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where £o is the coherence length, iV(0) is the density 
of normal electron states at the Fermi surface, /' is the 
slope of the order parameter at the vortex axis r — 0, I is 
a numerical constant, <I>o is the flux quantum and B(H) 
is the magnetic induction. The factor B(H)/&q is the 
density of vortices. The coherence length in Eq. (0) is 
related to the value of the energy gap Ao far outside the 
vortex cores i.e., £o = Jt-vf/tvAq. In Ginzburg-Landau 
(GL) theory, / w r/f at r = 0, so that /' « l/£ 
The GL coherence length £ is the characteristic length 
scale for spatial variations in the order parameter, and 
unlike £o> depends on temperature and is proportional 
to the field-dependent vortex-core radius. Using Eq. (Q), 
the contribution of the vortex cores to the specific heat 
when fc B T> ^l/E F is 



S (T,H) = -7v 2 N(E)k 2 B T = Tv 2 £ 0ln T 



B(H) 



(2) 

where 7„T = (2 / i)iv 2 N (Q)k 2 B T i s the normal-state elec- 
tronic specific heat and Ep is the Fermi energy. Thus, 



the specific heat in the vortex state may be written as 

'B(H) 



C(T,H) = C(T,0)+7v^olnT 



(3) 



where C(T, 0) is the specific heat at zero field in the 
Meissner state — with values from Fig. 1 given in Table @. 

Maki p(J] was the first to consider the effect of vortex- 
vortex interactions on the specific heat — deriving the fol- 
lowing relation from London theory, just above H c \ in the 
high-K limit 

TBd 2 H cl , TB 3 / 2 \ UH cl x 1 
H - H cl \ dT 

1 d\ 




XdT ' 



(4) 



where A is the magnetic penetration depth. Strictly 
speaking, B and A are functions of both T and H . The 
second term in Eq. (|]) has the same form as Eq. (Q), 
whereas the third term describes the perturbation due to 
vortex-vortex interactions. In Ref. jll|] it was suggested 
that this latter term may explain the downward curva- 
ture in C(H) observed at low fields in V3SL However, 
Eq. ^ has never been quantitatively verified. Further- 
more, it was noted in earlier measurements of C(T,H) 
in V 3 Si (Ref. [|l|), that Eq. (|) does not yield a term 
~ HT 3 which is observed experimentally. 

Recently, Ichioka et al. used the quasiclassical 
Eilcnberger equations to calculate the DOS in an s-wave 
superconductor with interacting vortices. At T = K the 
DOS is found to vary as H 67 due to a change in the 
slope of the order parameter in the vortex cores. Fitting 
the increasing field ZFC data for C(T, H)/T in the inset 
of Figs. 3 and 4 (the fit is not shown) to a function of the 
form C(T, 0)/T+(3H n using the values of C(T, 0)/T from 
Fig. 1, gives (3 = 5.36(4) mJ/mol K 2 T™ and n = 0.75(1) 
at T = 2.3 K, and /3 = 4.70(2) mJ/mol K 2 T" and 
n = 0.826(8) at T = 4.2 K. Using these two values 
of the exponent n, a linear extrapolation to T = K 
gives n = 0.66 — which is close to the result predicted in 
Ref. p2[ . Thus within this theoretical picture, vortex- 
vortex interactions affect the specific heat by decreasing 
both the size of the vortex cores and the corresponding 
density of bound QP states. 

A direct test of this idea can be made using the pre- 
cise field dependence of the vortex-core radius deter- 
mined from /LtSR. We first note that for a fixed value 
of the temperature, the radius of a coherence length £ 
in Eq. (||) is usually assumed to be independent of H, so 
that £ ~ £0 — in which case the relation $0 = 2tv£ 2 H C 2 can 
be used to give the familiar result: C corC s ~ JnTB / H C 2- 
More generally, however, the density of core states N(E) 
in Eq. (Q) will have two sources of magnetic field de- 
pendence. First, the vortex density B(H)/$q increases 
rapidly just above H c i approaching H/$q as H H C 2- 
Second, due to vortex-vortex interactions, ro ~ £ de- 
creases rapidly just above H c \ approaching a constant 
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value near H c2 . In Ref. |T(|, the fitted value of the 
GL coherence length for NbSe2 is well approximated by 
the relation £(H) = 46(2) + 28.9(9)/(ff - # c i) 1/2 and 
£(if) = 47(5)+46(3)/( J H-iT c i) 1 / 2 at T = 2.3 K and 4.2 K, 
respectively. Substituting these relations into Eq. (j^) 
gives an equation of the form 

C(T,H) = C(T,0)+ Cl TB(H) + T^f^ , ( 5 ) 

where ci and c 2 are numerical constants. The form of 
Eq. (||) is similar to Maki's macroscopic specific heat 
equation, since d 2 H c i/dT 2 < in Eq. (Eh. The main 
difference is that the second term in Eq. ^) was derived 
from data at fields 4JJ c i < H < 0.3H c2 , whereas Eq. (Q) 
is valid near H c i only. 

The form of B(H) in NbSe2 was determined from the 
ZFC magnetization measurements with increasing mag- 
netic field, where B(H)—H+AttM. A phenomenological 
equation was chosen to fit the magnetization data (e.g., 
solid curve in Fig. 2, inset) since there is no analytical 
equation from theory which is valid over the entire field 
range. We find that 

B(H) « J ff-0.00009(2) J ff" 1 - 32(7) 

- 0.0062(2) ln(0.91(6)/if) , 
B(H) w iJ-0.00063(2)ff-° 65(7) 

- 0.00105(12) ln(0.87(8)/#) 

at T = 2.3 K and 4.2 K, respectively. The contribution 
of B(H) to the field dependence of C(H) is observed 
by assuming £ is independent of field and equal to f;o> 
and substituting the relation for B(H) into Eq. (||). The 
solid curve in the inset of Figs. 3 and 4, is the corre- 
sponding fit which gives C(T,0)/T = 9.2(1) mJ/mol K 2 
and 30.25(5) mJ/mol K 2 at T = 2.3 K and 4.2 K, re- 
spectively. It is clear that the field dependence of the 
magnetic induction above H c \ does not alone account 
for the low- field curvature in C(H). We note that the 
fitted value of C(T,0)/T at T = 2.3 K, is nearly 20 % 
larger than that measured in zero field. 

To account for the effect of vortex-vortex interactions, 
we substitute the precise expression for £(T, H) measured 
by ^SR fjj| into Eq. (|) to give 

C(T,H) C(T,0) m [l + &i(r) H] 

—f— = —F~ + a{T)B{T > H) [1 + b 2 (T) H] ' (6) 

where a(T) = 7r 2 £ 7n£(r, 0)//$ , h(T) and b 2 (T) are 
temperature dependent coefficients. The values for b\(T) 
and b 2 (T) were determined in the work of Ref. Jill and 
are given in Table ffl. Note that £(T, 0) is merely an ex- 
trapolation of £(T, H) from the vortex state to H = T, 
and should not be confused with the coherence length in 
the Meissner state. The unusual form of the fiSK term 
has the following interpretation: (?) the numerator is the 
field dependence of the fitted magnetic penetration depth 



\ a b whereas, (ii) the denominator is the field dependence 
of the GL parameter K=\ a b/£ a b- 

In the main panel of Figs. 3 and 4, the data is plotted 
over the field range of the fj,SK experiment. The solid 
curve is a fit to Eq. (^), with the fitted values of a(T) 
and C(T, 0)/T given in Table Q. The quality of these fits 
lend strong support to the assertion that the downward 
curvature of C(H) at low fields is partially due to the 
expansion of the vortex cores. Using the fitted values 
of a(T), 7„ = 15.7(3) mJ/mol K 2 from Fig. 1, the GL 
value I = 1.92 §, £(2.3, 0) = 157(4) A and f (4.2, 0) = 
252(10) A from Ref. (Tefl, we find that £ = 56(4) A and 
Co = 42(8) A at T = 2.3 K and 4.2 K, respectively. 

As discussed earlier, in a pure d K 2_ y 2-wave supercon- 
ductor the DOS comes predominantly from the QP spec- 
trum outside of the vortex cores, since the presence of 
the gap nodes inhibit the formation of localized states 
within the cores. In this case the low-field expansion 
of the vortex cores recently measured by /iSR in YBCO 
p3| , p4j , should have little effect on C(H). However, the 
vortex-vortex interactions which are responsible for the 
changing core size should play some role in the behavior 
of C(H). For instance, in Ref. |22j it was shown that 
due to nearest neighbor vortex-vortex interactions, the 
density of the extended QP states in a d x 2_ y 2~wa.\e su- 
perconductor will be proportional to H° , rather than 
~ VH, as predicted by Volovik M. Thus far, measure- 
ments of C(H) in the HTSCs have not been performed 
with sufficient accuracy to verify this prediction. 

In conclusion, we have shown that the field dependence 
of the electronic specific heat measured in the s-wave su- 
perconductor NbSe2 is due to both the field dependence 
of the magnetic induction and the vortex-core size. We 
attribute the latter contribution to vortex- vortex interac- 
tions. Our findings provide a simple explanation for the 
downward curvature of C(H) reported in superconduc- 
tors which do not have nodes in the energy gap function. 

We thank A.V. Balatsky, L. Taillefer and M. Chiao 
for helpful and informative discussions. The work at Los 
Alamos was performed under the auspices of the US De- 
partment of Energy. 
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TABLES 



T 
[K] 


C(T,Q)/T 

r mj ] 


a 

r mj i 


6^ 


IT" 1 ] 


L mol K' J i 

Fig. 1 Eq. (pi) 


Lmol K^Tj 


2.3 
4.2 


7.7(1) 7.78(4) 
29.6(1) 29.26(9) 


12.7(1) 
15.0(5) 


0.555 
0.891 


2.48 
4.68 



TABLE I. Results from fitting the specific heat for NbSe2 
to Eq. (|). Only C(T, 0) and a(T) were free to vary. 



FIGURE CAPTIONS 

Figure 1. The specific heat of NbSe2 in zero field and 
H = 4 T (i.e. the normal state), plotted as C(T, 0)/T 
vs T 2 . The solid line is a fit to C n /T = ln + (3T 2 , where 
D = (l2ir 4 Nk B /5/3) 1/3 is the Debye temperature. 

Figure 2. The specific heat of NbSe 2 at T = 2.3 K, 
plotted as C(T,H)/T vs H. The data was taken four 
different ways, as described in the text. The inset shows 
the DC magnetization-hysteresis loop for the ZFC crys- 
tal. 

Figure 3. The specific heat of NbSe 2 at T — 2.3 K 
over the field range of the /xSR experiment Q plotted 
as C{T,H)/T vs H. The solid curve is a fit to Eq. (g). 
Inset: The specific heat over an extended field range. 
The solid curve is a fit to Eq. (B|) . 



Figure 4. 
T = 4.2 K. 



Same as Fig. 3, but with the sample at 
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